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U ■ ABSTRACT 

-4— » , 

C/3 . 

, Context. Since 1999, we have been conducting a radial velocity survey of 179 K giants using the Coude Auxiliary Telescope at UCO/Lick 
• observatory. At present ~20-100 measurements have been collected per star with a precision of 5 to 8 m s"'. Of the stars monitored, 145 (80%) 
I show radial velocity (RV) variations at a level >20 m s"', of which 43 exhibit significant periodicities. 
^ Aims. Our aim is to investigate possible mechanism(s) that cause these observed RV variations. We intend to test whether these variations are 
intrinsic in nature, or possibly induced by companions, or both. In addition, we aim to characterise the parameters of these companions. 
Methods. A relation between log g and the amplitude of the RV variations is investigated for all stars in the sample. Furthermore, the hypothesis 
that all periodic RV variations are caused by companions is investigated by comparing their inferred orbital statistics with the statistics of 
companions around main sequence F, G, and K dwarfs. 
I Results. A strong relation is found between the amplitude of the RV variations and log g in K giant stars, as suggested earlier by Hatzes & 
' Cochran (1998). However, most of the stars exhibiting periodic variations are located above this relation. These RV variations can be split in a 
' periodic component which is not correlated with log g and a random residual part which does correlate with log g. Compared to main-sequence 
, dwarf stars, K giants frequently exhibit periodic RV variations. Interpreting these RV variations as being caused by companions, the orbital 
1^ parameters are different from the companions orbiting dwarfs. 
. ^ Conclusions. Intrinsic mechanisms play an important role in producing RV variations in K giants stars, as suggested by their dependence 
on log g. However, it appears that periodic RV variations are additional to these intrinsic variations, consistent with them being caused by 
companions. If indeed the majority of the periodic RV variations in K giants is interpreted as due to substellar companions, then massive 
planets are significantly more common around K giants than around F, G, K main-sequence stars. 
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1. Introduction 

For more than a decade, radial velocity observations with ac- 
curaci es of order ms~ ' have bee n within reach (see for in- 
stance iMarcx^Butie] ( I2OOOI) an dlOueloz et all (l200lb ). Even 
accuracies of less than 1 ms~' dPepe et al.l l2003h are possi- 
ble now. With these observations, more than 200 sub-stellar 
companions have been discovered by measuring the reflex mo- 
tions of their parent stars. Most of these sub-stellar compan- 
ions have been detected around F, G and K main sequence 
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stars, but detections ar ound an A sta r ( Galland et alj 200^ 



(l2006h . 



Johnson et al 



and several subgiants dJohnson et al. 
(2()07)) have also been reported recently. Moreover, 10 
giant stars were re ported to have su b-stellar companions 
(t Dra conis (K2III) |FnnkeLalJ (l2002h. HP 104985 (09111) 

(i2003h . 
HD 11977 (G5III ) 
20061), 
(2007), 



Draconis (K2III) iFrink et al 
Sato et al. (20^, H D47526 (Kl III) 



HP 1 3 1 89 (K2 I I-IID Hatzes et al 



d2005[ 



Setiawan et al 



n 



ISetiawanet al.' ('2005'), Pollux (KOIII) 'Hatzes et al 
Reff'ert et al. (2006), 4UMa (KlIII) polUnger et a l 
NGC2423 No3 and NGC4349 Nol27 [Lovis & Mayor! d2007l) 
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and recently HD 17092 (KOIII) iNiedzielski et all (l2007b fl. In 
addition to searches for extra-solar companions, radial veloc- 
ity observations prove to be very useful for detecting solar- 
like oscillations in stars with turbu lent at mospheres, such as the 
dwa rf a Cen A (e.g. Adding et al.ll2006l). t he su bgiant Procyon 
(e.g.'Eggenberger et al. 2004: Martic et all2004i) and the giant 
€ Ophiuchi (e.g. De Ridder et al. 200.^^ 

With techniques for accurate radial velocity observations 
at hand, a survey was started in 1999 to verify whether K 
giants are stable enough to be used as astrometric refer- 
ence stars for SIM /PIanetQuest (Space Interferometry Mission) 
( Frink et aDl2001 ). This survey contains 179 stars and uses the 
Coude Auxiliary Telescope (CAT) at University of California 
Observatories / Lick Observatory, in conjunction with the 
Hamilton Echelle Spectrograph. The survey has recently been 
expanded to about 380 giants and is still ongoing. For the anal- 
ysis described in the present paper only data from the initial 
179 stars are used. 

From this survey, comp anions have b een announced for t 



2006h . 



Draconis (IFrink et alJ 120021) and Pollux (IReffert et al. 
Stars with radial velocity variations of less than 20 ms"' have 
been presented as stable stars bv Hekker et al. (2006). In addi- 
tion, some binaries discovered with this survey, as well as an 
extensive overview of the sample, will be presented in forth- 
coming papers. 

As almost all of the stars show significant radial velocity 
variations, we investigate here which mechanism causes these 
variations. Non-periodic radial velocity variations, of the order 
of the investigated timescales, are most likely caused by some 
intrinsic mechanism, while the periodic variability can also be 
caused by companions. We also investigate the characteristics 
of these companions. 

In Sect. 2, the radial velocity observations are described in 
detail. In Sect. 3, the relation between the observed radial ve- 
locity amplitude and surface gravity is investigated. In Sect. 4, 
we explore the hypothesis that all periodic radial velocity vari- 
ations are caused by sub-stellar companions, and we compare 
the inferred orbital parameters with those obtained for sub- 
stellar companions orbiting main sequence stars. Our conclu- 
sions are presented in Sect. 5. 

2. Radial velocity variations 

The initial 179 stars selected for the radial velocity survey are 
used in the present work . These stars h ave been selected from 
the Hipparcos catalo gue dPerrvnian & ESA.1997) . based on the 
criteria described by Frink et al. ( 200 lb . The selected stars are 
all brighter than 6 mag, are presumably single and have pho- 
tometric variations < 0.06 mag in V. The survey started in 
1999 at Lick Observatory using the Coude Auxiliary Telescope 
(CAT) in conjunction with the Hamilton Echelle Spectrograph 
(R=60000). The syst em with an iodine cell in t he light path 



has be en described bv lMarcv & Butlen (Il992h and lValenti et al 



( 119951) . With integration times of up to thirty minutes for the 
faintest stars (nii, = 6 mag) we reach a signal to noise ratio of 
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Fig. 1. Radial velocity variations as a function of phase for a 
star (HIP34693) with a highly significant period (top), with its 
periodogram (bottom). The dashed line in the periodogram in- 
dicates the significance threshold. 



about 80 - 100 at /I = 5500 A, yielding a radial velocity pre- 
cision of 5 - 8 ms"'. As we are looking for radial velocity 
variations of order 10 to 100 m s"', this is adequate and hence 
no attempt has been made to reach the 3 m s~' accuracy w hich 
is in principle possible with this setup (iButler et al. 1996). For 
the determination of th e radial velocities the pipeline described 



by iButler et al.l (11996 ) is used. In this pipeline, a template io- 
dine spectrum and a template spectrum of the target star ob- 
tained without an iodine cell in the lightpath are used to model 
the stellar observations with a superposed iodine spectrum. The 
Doppler shift is a free parameter in this model. Note that with 
this method the absolute radial velocity is not measured, but the 
radial velocity relative to the stellar template is obtained. 

All 179 stars are subjected to a period search. The peri- 
odicity of the radial velocity variations is determined first of 



all fro m a classical Lomb-Scargle (LS) periodogram JScargk 



' For updated information on sub-stellar 
|http://exoplanet.eu and littp ://exo planets . org [ 



companions, see 



19821) . The significance threshold is set to 6cr, where the noise 
level is determined from the average power of the residual 
Scargle periodogram for frequencies between and 0.03 cy- 
cles per day (c/d) (0.35 //Hz) and a frequency step of 0.00001 
c/d (0.12 • 10"^ //Hz). We adopted the conventional method of 
iterative sinewave fitting ('prewhitening ') to search for subse- 
quent frequencies d Kuschnig et al. 1997 ). In Figs [T] and |2] the 
radial velocity variation as a function of phase is shown for 
two stars. The period of the star in Fig.[T]is highly significant, 
while the one in Fig. |2]is close to the significance threshold. 
Periodograms are shown in the bottom panels of these figures. 
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Table 1. Single stars with significant frequencies. 
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Fig. 2. Radial velocity variations as a function of phase for a 
star (HIP7607) with a period close to the significance threshold 
(top), with its periodogram (bottom). The dashed line in the 
periodogram indicates the significance threshold. 



As properly emphasized bv lCumming et al. (Il999l) . such a clas- 
sical period search may not be appropriate for unevenly spaced 
sparse data, even though we set the significance level at a con- 
servatively high level. In order to check this, we have done an 
additional LS analysis after prewhitening the original data by a 
linear polynomial. This led almost always to the same frequen- 
cies. We only accepted a frequency when it was found to meet 
the significance criterion for both these analyses. The signifi- 
cant frequencies are listed in Table[T] 



3. Radial velocity amplitude 
relation 



surface gravity 



Hatzes & CochranI ( 119981) already investigated the origin of the 
observed radial velocities in K giant stars. Although their sam- 
ple contained only 9 stars, they suggested that the amplitude 
of the radial velocity increases with decreasing surface gravity 
(logg). In lower surface gravity it takes longer to decrease the 
velocity of a moving parcel which results in larger amp litudes 
and the relation suggested by lHatzes & CochranI (1 1998h would 
therefore be evidence for pulsations or rotational modulation as 
the mechanism for these long period radial velocity variations. 

For the pres ent sample, log g values were determined spec- 
troscopically by iHekker & MelendezI ( l2007h . by imposing ex- 
citation and ionisation equilibrium of iron Unes through stel- 
lar models. The equivalent width of about two dozen care- 
fully selected iron lines were used for a spec troscopic LTE 
analysis based on the 2002 version of MOOG (ISnedenlll973h 
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and Kurucz model atmospheres which include overshooting 
( Castelli et al.ll997 ). These authors estimated the error on log g 
to be 0.22 dex from the scatter found in a comparison with lit- 
erature values. A detailed description of the stellar parameters 
for individu al stars and a compar i son w ith literature values are 
available in iHekker & MelendezI (l2007b and is therefore omit- 
ted here. 

In Fig. [3] we show half of the peak-to-peak value of the 
observed radial velocity variations as a function of log g for K 
giants in our sample. A clear trend is visible between increasing 
radial velocity variations in single stars, and decreasing \ogg, 
which provides a strong indication that, at least for a large frac- 
tion of stars in our sample, the observed radial velocity vari- 
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Fig. 3. Half of the peak-to-peak variation of the radial velocity as a function of surface gravity (logg). The blue o indicate the 
single stars with periodic radial velocity variations (see text for periodicity criteria), the red □ indicate single stars with random 
radial velocity variations, and stellar binaries (companion mass > 100 Mjup) are indicated with black o symbols. (Colours are 
only visible in the online version of the paper.) The solid line is the best fit through the random single stars, the dotted line 
indicates the Icr interval around the best fit and the dashed line indicates the 3.5cr interval. Single stars with a higher radial 
velocity amplitude than expected based on their logg value (more than 3.5cr above the best fit) are indicated by arrows. Six of 
the 8 stars with periodic radial velocity variations and logg < 1.6 are classified bright giants or supergiants ( Ferryman & ES aI 



1997). 



ations are induced by a mechanism intrinsic to the star. This 
trend is present for stars with random as well as stars with peri- 
odic radial velocity variations. Also, nearly all stars with peri- 
odic radial velocity variations and log g > 1 .6 are located above 
the fit in Fig. [3] Seven single stars have a higher radial veloc- 
ity variation than expected based on their logg value, i.e. they 
are situated more than 3.5cr above the best fit for the relation 
obtained for single stars. These stars are indicated with arrows 
in Fig. [3] The radial velocity variation observed for H1P53229 
may be due to a stellar companion in a wide orbit, with a period 
much longer than the observation time span. Due to this long 
period the companion mass, and, therefore, the (sub-)stellar na- 
ture, is still very uncertain. H1P33152 is classified as a super- 
giant. The observed radial velocity variations for H1P80693, 
H1P36616, and HIP88048 can be fitted with two Keplerian or- 
bits, while HIP75 458 can be explai ned by an eccentric sub- 
stellar companion dFrink et al J 120021) and an additional linear 
trend, indicating a companion in a wide orbit. H1P34693 can be 



fitted very accurately with a Keplerian orbit of a single nearly 
sinusoidal sub-stellar companion. 

In order to investigate the simultaneous occurrence of sub- 
stellar companions and oscillations in giants, the best Keplerian 
fits are subtracted from the observed radial velocity variations 
of both the binaries and the single stars with significant periodic 
radial velocity variations. The half peak-to-peak values of these 
residuals are plotted as a function of surface gravity in Fig.H] 
also showing the linear fit through the non-periodic stars, and 
the amplitudes of the subtracted periodic signals. There are sev- 
eral interesting points to make about this graph. Firstly, almost 
all periodic stars show larger radial velocity variations than pre- 
dicted by the relation found for non-periodic stars, but when the 
periodicities are removed, their residual radial velocity varia- 
tions follow the same relation as found for the non-periodic 
stars. This could be interpreted as evidence for both intrinsic 
(non-periodic radial velocity variation) and extrinsic (periodic 
variations) mechanisms playing a role in these stars. Secondly, 
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Fig. 4. Top: Half of the peak-to-peak variation of the radial velocity as a function of surface gravity (logg), as in Fig. [3] but 
showing only those stars with periodic radial velocity variations (blue o) and stellar binaries (black o). The blue + and black x 
symbols indicate the amplitude of the radial velocity variations for these stars after subtraction of the Keplerian fits. (Colours are 
only visible in the online version of the paper) The solid line indicates the linear fit through the stars with non-periodic radial 
velocity variations (from Fig. [3]). Bottom: Amplitudes of the subtracted Keplerian fits as a function of logg. 



there is no correlation between the amplitude of the subtracted 
periodic signal and log g, which provides additional evidence 
for the presence of companions. Thirdly, almost all (8 out of 9) 
stars with logg < 1.6 exhibit periodic variations. If these in- 
deed have an extrinsic mechanism, it would mean that ~ 90% 
of these stars have sub-stellar companions. However, from an 
astrophysical point of view, stars with such low surface grav- 
ities are already very high on the giant branch or even on the 
asymptotic giant branch. At these low gravities, stars cannot be 
constant anymore, as the outer layers are so diluted that insta- 
bilities occur easily, either periodic or random. Stated differ- 
ently, these stars are very close to the semi-regulars, which are 
on their way to become Mira variables. Hence, these periodic 
variations could well be intrinsic. 



4. Companion Interpretation 

From the analysis of the correlation between radial velocity 
amplitude and surface gravity we have evidence for the pres- 
ence of both intrinsic variability and companions in at least a 
fraction of the K giant stars with periodic radial velocity varia- 
tions. In order to study the characteristics of these companions, 
we take the hypothesis that the periodic radial velocity varia- 
tions detected in 43 of our K giant stars, excluding binaries, are 
caused by sub-stellar companions. Under this hypothesis, we 
investigate the statistical properties of the orbital parameters of 
the sample and compare these with the statistical properties of 
companions orbiting F, G and K dwarfs. 

According to our analysis, 55 stars in the sample would 
have a single companion, and 11 stars multiple companions. 
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Twenty-three (22 single and 1 in a multiple system) of these 
companions have msini larger than 100 Mjup and should be 
interpreted as stellar binaries. By advancing the multiple sub- 
stellar systems forward in time via a Runga-Kutta integration, 
we investigated the stability of the systems, taking into ac- 
count the mutual interactions of the companions. With the or- 
bital parameters that minimise taken at face value, we found 
that most of the inferred sub-stellar multiple systems would be 
'likely unstable', with a change in semi-major axis > 1% and 
< 10%, or 'unstable' with a change in semi-major axis > 10% 
on a time scale of 100 years due to companion-companion in- 
teraction. 

However, the inferred stability depends on the starting 
epoch of the computations, as well as on the orbital parame- 
ters, which might change with an increasing number of obser- 
vations. Furthermore, there is no guarantee that the obtained 

minimum is a global minimum. Therefore, stars with mul- 
tiple inferred sub-stellar companions that seem to be unsta- 
ble, might also have stable solutions. One c ould also use the 
equ ations for dynamical s tability de s cribed bv lGladmanI j 19931) 
and iMarcv et all (l200lh . iGladmanI (Il993h also notes that the 
Hill stability criteria for companions in initially eccentric or- 
bits may not be met, but that the systems may still be found to 
be empirically quite stable for a long period of time. In order to 
draw a firm conclusion on the stability of a particular system, 
a more thorough investigation is needed, as well as data with a 
longer time base, which is beyond the scope of this paper 

For all stars with p eriodic radial velocity v ariations, we 
checked the Hipparcos (IPerryman & ESAlll997n photometry. 
We checked periodograms for significant frequencies close to 
the obtained radial velocity period, and plotted the photomet- 
ric values phased with the radial velocity period. None of the 
stars show photometric variations related to the observed radial 
velocity variations. 

The mass distribution of our K giant sample is not known 
very well. The stellar masses are typically between 1 and 4 
Mq. Hence most of their main sequence progenitors should 
have been of A or F spectral class. The distribution of orbital 
parameters of sub-stellar companions orbiting A and F main 
sequence stars is still unknown. The core accretion model pre- 
dicts more giant planets around more massive stars, so that the 
distribution of orbital parameters of sub-stellar companions or- 
biting F, G and K main-sequence stars probably cannot serve as 
a proxy. However, it should be instructive to compare the two 
distributions. 

Since the data presented here span ~ 2500 days, radial ve- 
locity variations with longer periods are uncertain, and, there- 
fore, not taken into consideration. Companions with periods 
exceeding the observation time span are also excluded from 
the F, G, and K main sequence star statistics. 

4.1. Mass distribution 

Fig. |5] shows the distribution of inferred companion masses 
of our K giant sample. First, notice that 30% of the inferred 
companions would have masses in the brown dwarf regime 15 
Mjup < m sin i < 80 Mjup. This is in sharp contrast to the brown 



Companion Assumption 




msini \M , 1 



Fig. 5. A histogram of m sin / of inferred companion masses or- 
biting K giants in our sample. The dashed line indicates the 
rise of planet masses M"' "^ from 10 Mjup down to Saturn 
masses for sub-stel lar companions around main sequence stars 
(IMarcv et al.ll2005h . normalised to the number of stars in our 
sample. 
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Fig. 6. Zoom in on the low end of the companion mass distri- 
bution of inferred sub-stellar companions orbiting K giant stars 
in our survey, shown as percentage of the total number of stars 
in the sample (gray histogram). The hatched histogram shows 
the distribution of compa nion masses orbiting main sequence 
stars as shown in Fig. 1 of lMarcv et al.l (1200 5). as a percentage 
of the total number of stars in their sample. The cross hatched 
area are main sequence stars with companions at a semi-major 
axis smaller than 0.3 AU (see Fig.|7]i. 



dwarf statistics around F, G and K main sequence stars for 
which only very few companions are found with m sin / > 15 
Mjup around more than thousand stars. This is known as the 
brown dwarf desert and is possibly caused by migration and 
merging of brown dwarfs in a viscous dis k with a mass at least 
comparable to the brown dwarf mass ( Armitage & Bonriell 
2002). 

The dashed line in Fig. |5] indicates the rise of sub-stellar 
companion masses M""'^ fr om 10 Mjup dow n to Saturn 
masses for main sequence stars dMarcv et al 12005 ), normalised 
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Fig. 7. Semi-major axis distribution of inferred K giant com- 
panions in our sample (gray histogram), shown as a percentage 
of the total number of stars in the sample. The hatched his- 
togram shows the distributi on for F, G and K main sequence 
stars as shown in Fig. 2 of iMarcy et al. ( 20051) . as a percent- 
age of the total number of stars in their sample. The main se- 
quence stars with a companion orbiting at a semi-major axis 
smaller than 0.3 AU are cross-hatched. These are also indicated 
inFigs|6]and[8l 



to the number of stars in our sample . We use a two-sided 
Kolmogorov-Smirnov test dPress et al. 1992 ) (hereafter KS- 
test) to compare the M"'°^ fit and the mass distribution 
(m sin / < 80 Mjup) of our sample and find a probability of 
0.05% (0.002% for msini < 28 Mjup) that these are identi- 
cal. This implies that inferred sub-stellar companions around 
K giants in our sample have higher masses compared to com- 
panions around F, G and K main sequence stars. 

Fig. |6] shows the low-mass companion distribution of our 
survey. Most K giant companions would have inferred masses 
between 2 and 8 Mjup, while the fraction of companions or- 
biting F, G and K dwarfs strongly decreases with increasing 



4.2. Semi-major axis distribution 

The distribution of the companions' semi-major axis is shown 
in Fig.|7] No inferred companions with semi-major axis smaller 
than 0.3 AU are present in the K giant sample, possibly due to 
increased stellar radii of giants. The fraction of stars with an in- 
ferred companion with a semi-major axis between 1 and 3 AU 
is much higher among the K giants compared to the F, G and K 
dwarfs. A comparison between the two distributions with a KS- 
test reveals a probability for the two distributions to be iden- 
tical of 11%. This increases to 32% when omitting the main 
sequence stars with semi-major axis < 0.3 AU. The increasing 
incompleteness beyond 3 AU, due to the limited time span of 
surveys, is present in both samples as the surveys cover a com- 
parable amount of time. This incompleteness cannot cause the 
significant diff'erence in the peak between 1 and 3 AU. 



Fig. 8. Period distribution of the observed radial velocity vari- 
ations shown as a fraction of all significant periods (gray 
histogram). The period distribution of main sequence stars 
( Butler et al..2006..) is shown in the hatched histogram, as a frac- 
tion of the total number of observed companions. The cross 
hatched area are main sequence stars with a semi-major axis 
smaller than 0.3 AU (see Fig.|2l). 



4.3. Period distribution 

In Fig. |8] the period distribution of the observed radial veloc- 
ity variations is shown and compared with the companion pe- 
riod distribution of dwarfs. The large fraction of F, G and K 
dwarf companions with orbital periods shorter than 100 days 
corresponds to the ones with semi-major axis smaller than ap- 
proximately 0.3 AU. The close-in short-period companions are 
not present around K giants, while about 80% of these stars 
with observed radial velocities have periods ranging between 
400 and 800 days. A KS-test reveals a probability of less than 
0.0001% for the two distributions to be identical. The proba- 
bility remains below this level, when the companions orbiting 
main sequence stars with a semi-major axis < 0.3 AU are omit- 
ted. 

4.4. Eccentricity distribution 

Fig. |9] shows the distribution of companion eccentricities for 
K giants in our sample and for dwarfs. Companions of dwarfs 
with periods less than 20 days are excluded, as these might be 
tidally circularised. The fraction of companion eccentricities 
< 0.3 for the giants is 75% compared to 50% for companions 
orbiting F, G and K dwarfs. The KS-test shows that these dis- 
tributions are nearly identical (97%). 

4.5. Iron abundance 



Companion occurrence correlates strongl y with the abun - 
dance of heavy e lements (se e for instance Gonzalea (Il997h . 
Fischer & Valentil §005) and ISantos et all ( l2005h ). such that 
F, G, and K dwarf stars with supersolar abundance are more 
likely to harbour sub-stellar companions (about 50% of the 
stars with 0.3 < [Fe/H] < 0.5). The increase of the fraction 
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Fig. 9. Distribution of eccentricities for possible companions 
around K giants in our sample (gray histogram) shown as a 
fraction of all possible sub-stellar companions in the sample. 
The hatched histogram is the eccent ricity distribution o f com- 
panions around main sequence stars (iButler et al ]l2006h shown 
as a fraction of all companions around main sequence stars. 
Companions with periods shorter than 20 days are excluded 
from the latter sample as these might be tidally circularised. 



of F, G and K dwarfs harbouring companions with increasing 
metallicity is well fitted with a power law, yielding a prob- 
ability for such a star to harb our a companion to b e: P = 



0.03 ■ [(A^Fe/A^H)/(A^Fe/A^H)0]^ " ( Fischer & Valentill2005 



In Fig. [To] the iron abundance distribution of stars with pe- 
riodic radial velocity variations is shown as a percentage of the 
total number of observed stars with iron abundance in the same 
interval. The iron abundance is determined spectroscopically 
by imposing excita tion and ionisation equilibri um in iron lines 
and is described by lHekker & Melendea (120071) . The maximum 
iron abundance of a K giant star in our sample is 0.29 and, 
therefore, we do not probe the high metallicity region in which 
F, G and K dwarfs are most likely to harbour a companion. 

The mean metallicity of the K giants in the entire sample is 
-0.12 dex, while the mean metallicity of the stars with periodic 
radial velocity variations, presented in Fig. [TO] is -0.13 dex. 
No correlation between companion occurrence and abundance, 
similar to the one which is present for dwarf stars, is found in 
this sample of giant stars. 

5. Discussion and Conclusion 

The tight correlation we found between log g and half of the 
peak-to-peak radial velocity variations seems to indicate that a 
large fraction of the observed radial velocity variations in our 
sample of K giants is induced by mechanism(s) intrinsic to the 
stars. We also present evidence that both intrinsic and extrinsic 
mechanisms play a role. The stars with a significant periodic 
signal are almost all located above the radial velocity ampli- 
tude vs. log g relation, but when the periodic signal is removed, 
the residuals show the same trend as for the non-periodic stars. 
Furthermore, no correlation is present between the amplitude 
of the periodic signal and logg. 



Fig. 10. Iron abundance [Fe/H] distribution of K giant stars 
with periodic radial velocity variations shown as a percentage 
of the total number of observed stars with iron abundance in 
the same interval. The numbers above each bar on the his- 
togram indicate the ratio of stars with a significant periodic 
radial velocity variation to the total number of stars in each 
bin. The error bars are calculated assuming Poisson statistics 
(i.e., the percentage of stars with periodic radial velocities di- 
vided by the square root of the number of stars with periodic 
radial velocities). The dashed line is the power law derived for 
the increasing trend in the fraction of stars with companions 
as a function of metallicity of F, G and K main sequence stars 
dFischer & Valenti.2005,) . 



Almost all of the lowest logg stars show periodic varia- 
tions. It may be possible that stars with such low surface gravity 
cannot be constant and that in these dilute atmospheres insta- 
bilities can occur very easily, and therefore may be periodic, 
but not extrinsic. 

Based on the evidence that extrinsic mechanism(s) play a 
role for K giant stars with periodic radial velocity variations we 
investigated the hypothesis that this periodic signal is caused 
by the reflex motion of sub-stellar companions orbiting these 
stars. We presented the characteristics of the orbital parameters 
of these companions and compared them with the known or- 
bital parameters of sub-stellar companions orbiting F, G and K 
dwarfs. 

About 25% of the stars in our sample have radial veloc- 
ity variations with significant periodicity, and could possibly 
harbour a sub-stellar companion, while approximately only 
8% o f the 1330 F, G and K main sequence stars investigated 
by Marcv et al. ( 2005h ha ve a sub-stellar compan ion. Recently 
Ijohnson et al.1 (l2007h and iLovis & Mayod (l2007l) showed that 
the number of companion harbouring stars increases with mass. 
The giants in the present sample have typical masses between 
1 and 4 Mq and are in general more massive than the main se- 
quence stars investigated for companions. So, the high percent- 
age is qualitatively in agreement with the re sults from the lit- 
erature. Furthermore. lLovis & May on (120071) suggest that more 
massive stars form more massive planetary systems than lower 
mass stars. Figs |5] and |6] show that we find, in general, more 
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massive companions around the more massive K giants than 
are present around F, G and K dwarfs. 

The high percentage of more massive companions around 
the more massive K giant stars would also be compatible with 
the core accretion model. This model predicts very few giant 
planets, but a relatively large number of planet s with the mass 
of Neptune or sm aller around M dwarfs (Laughlin et al.ll2004 : 
Ida & Liiil 2005 ). This is mainly the result of a much reduced 
surface density of the disk and the resulting shorter disk evolu- 
tion timescales compared to those for more massive stars, im- 
plying that planet properties vary with the mass of their host 
stars. In particular, one would expect more sub-stellar compan- 
ions with higher masses in our giant sample, as is indeed the 
case if we assume that the companion hypothesis is correct. 

The mean metallicity of companion-hosting K giants 
would be similar to the mean metallicity of the total sam- 
ple. This is in contrast with the correlation between com- 
panion occu rrence and metallicity present in F, G and K 
dwarfs (e.g. iFischer & Valentil 120051) . So far, several groups 
have investigated the correlation between companion oc- 
currence and metallic ity f or giants with differ ent results. 
Sadakane et al.l (l2005b and IPasquini et al.l (l2007h agree that 
comp anion-ahosting giants are o n average not metal-rich, 
while Hekker & Melende3 (2007) find that giants with an- 
nounced companions have higher metallicities than their to- 
tal sample of giants. A det ailed discussion ab out these differ- 



ent results is presented by Hekker & Melendea (l2007h . They 
conclude that the samples on which the results are based 
are slightly different and the more m etal-rich stars used i n 
their study are lacking in the study by IPasquini et al.l (120071) . 
Furthermore, there is a difference in zero-point correction for 
the metallicities of announced companion-hosting stars from 
different surveys. All in all, these inferences are based on small- 
number statistics and all results have to be taken with caution. 

The larger semi-major axis and long periods of the inferred 
companions orbiting K giant stars compared to companions or- 
biting dwarf stars are most likely due to the extended atmo- 
spheres of K giants. For the eccentricity no significant differ- 
ence is found in the distribution between companions orbit- 
ing dwarfs or giants. Nevertheless, the high number of inferred 
companions around giants with eccentricities < 0.3 is striking. 
One could suspect that companion orbits circularise over time 
and that the companions in circular orbits are older than the 
eccentric ones, but there is no evidence for this hypothesis. 

In principle, nearly sinusoidal radial velocity variations 
could also be caused by pulsations or spots. Although the peri- 
ods of the radial velocity variations could well be the rotational 
periods of the stars, the presence of prominent spots is not very 
likely. In that case one would also expect photometric varia- 
tions with periods correlated wit h the radial velocity varia tions. 
From the Hipparcos photometry (IPerryman & ES All 19971) such 
correlations were not found. 

In order to distinguish with certainty between companions 
and pulsations as the cause of the observed radial velocity vari- 
ations, one needs to perform a spectral line profile analysis. 
A technique for doing this with very high-resolution spectra 
(R > 100000) will be presented separately (Hekker et al., in 



preparation) because the amount of such data at hand today is 
insufficient to add significantly to the conclusions of this paper 
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